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Linear  Feedback  Control  and  Numerical  Approximation 
for  a  System  Governed  by  the  Two-Dimensional  Burgers’  Equation 

R.  Chris  Camphouse,  James  H.  Myatt 


Abstract 

In  this  paper,  we  consider  the  problem  of  controlling  a  system  governed  by  a  two-dimensional  nonlinear  partial  differential 
equation.  Motivation  for  the  problem  is  the  development  of  control  methodologies  for  fluid  flow,  where  the  dynamics  of  the 
system  are  governed  by  the  nonlinear  Navier-Stokes  equations.  A  control  problem  for  the  two-dimensional  Burgers  equation 
is  considered.  An  initial  boundary  value  problem  is  formulated  to  model  a  right-traveling  shock  over  an  obstacle.  We  focus 
on  implementing  a  feedback  control  via  Dirichlet  boundary  control  on  the  obstacle.  Numerical  approximations  are  developed 
to  numerically  simulate  solutions  of  the  problem  with  and  without  control.  Numerical  examples  are  presented  to  illustrate  the 
efficacy,  as  well  as  the  shortcomings,  of  the  control  method.  Additionally,  the  influence  of  boundary  condition  on  the  functional 
gains,  and  the  resulting  controls,  will  be  demonstrated  through  numerical  examples.  Avenues  of  future  work  are  presented. 

Introduction 

One  of  the  commonly  studied  problems  of  flow  control  is  the  manipulation  of  flow  over  an  obstacle.  Typical 
control  objectives  include  velocity  tracking,  drag  reduction,  lift  increase,  and  separation  postponement  [13, 
16].  Due  to  the  complexity  of  the  Navier-Stokes  equations,  control  laws  are  often  determined  on  a  system-by¬ 
system  basis.  Experimental  data  are  taken  for  a  particular  system,  and  a  model  of  the  system  is  developed 
from  those  data.  Control  methods  are  then  developed  for  the  data-induced  model.  A  major  drawback  of 
this  approach  is  that  control  methods  developed  for  one  particular  conflguration  and  flow  condition  may  not 
directly  apply  to  a  different  configuration  and  flow  condition.  In  addition,  the  effectiveness  of  the  control 
developed  in  this  way  depends  on  the  accuracy  and  reliability  of  the  experimental  data.  By  their  very  nature, 
experimental  data  can  only  provide  a  sampling  of  the  flow  behavior  for  the  sampled  conditions.  As  a  result, 
dynamics  that  the  control  needs  to  suppress  or  enhance  may  be  missed  in  the  data  collection  process.  Ihe 
result  is  likely  to  be  a  control  that  is  not  as  effective  as  was  desired,  as  well  as  a  control  method  that  doesn’t 
generalize  to  different  conditions. 

Developing  control  methods  directly  from  a  partial  differential  equation  model  describing  flow  behavior 
has  many  advantages.  It  is  widely  accepted  that  the  Navier-Stokes  equations  provide  a  good  description  of 
fluid  behavior,  even  at  high  Reynolds  number.  As  a  result,  control  techniques  developed  from  the  Navier- 
Stokes  equations  have  the  possibility  of  being  relevant  to  a  variety  of  flow  conditions  and  configurations.  In 
addition,  general  optimality  results  can  often  be  obtained  when  formulating  controls  directly  from  governing 
equations,  rather  than  equations  obtained  via  sets  of  data. 

There  are  many  difficulties  inherent  in  developing  feedback  control  laws  for  the  Navier-Stokes  equations. 
One  difficulty  is  due  to  the  nonlinearities  involved.  There  is  an  extensive  number  of  control  techniques 
available  for  linear  distributed  parameter  systems.  However,  the  theory  needed  to  control  highly  nonlinear 
problems  is  far  from  complete.  Another  difficulty  is  that  numerical  simulation  of  fluid  flows  can  be  very 
expensive  due  to  the  large  number  of  states  needed  to  resolve  a  given  flow  field. 

The  scalar  Burgers’  equation,  which  has  a  nonlinearity  similar  to  that  in  Navier-Stokes,  is  often  used 
for  the  development  of  control  methods  relevant  to  flow  control  [7,  8,  11].  These  developments  are  often 
done  for  the  one-dimensional  case.  Optimal  actuator  and  sensor  placement,  as  well  as  the  control  objective, 
are  more  complicated  issues  to  consider  in  the  case  of  higher  dimensions.  In  this  paper,  we  implement  a 
feedback  control  for  a  system  whose  dynamics  are  governed  by  the  two-dimensional  Burgers  equation  over  a 
geometry  similar  to  those  of  interest  in  flow  control  problems  involving  an  obstacle.  Applying  these  control 
techniques  to  an  aerodynamic  flow  control  problem  is  the  subject  of  a  future  paper. 

Problem  Description 

In  this  section,  we  develop  a  partial  differential  equation  model  for  the  system  under  consideration.  To 
this  end,  let  C  he  the  open  rectangle  given  by  {a,b)  x  (c,d)  where  a,b,c,d  €  K.  Let  02  £  Le  the 
rectangle  given  by  [01,02]  ^  l^it  where  o  <  oi  <  02  <  h  and  c  <  bi  <  b^  <  d,  i.e.,  fl2  C  Oi-  The  problem 
domain,  O,  is  given  by  0  =  Oi  \i22-  In  this  configuration,  O2  is  the  obstacle  on  which  we  implement  Dirichlet 
boundary  control. 
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The  dynamics  of  the  system  are  given  by  the  two-dimensional  Burgers’  equation 


-u;(t,x,y)  +  Ki- 


(1) 


for  f  >  0  and  {x,  y)  €  Q.  Behaviors  described  by  Burgers’  equation  include  shock  formation,  shock  propa¬ 
gation,  and  rarefaction  wave  formation.  In  (1),  K\  and  K2  are  constants  used  to  scale  the  nonlinear  terms. 
The  quantity  Re  is  a  nonnegative  constant,  and  is  analogous  to  the  Reynolds  number  in  the  Navier-Stokes 
equations. 

In  order  to  fully  specify  the  model,  we  need  to  specify  conditions  on  ^17 1  and  ^^2,  as  well  as  an  initial 
condition.  As  we  are  interested  in  implementing  control  on  00.2^  denote  the  sides  of  dQ.2  as 


Ti  -  {(x,fei)  1  ai  <  x  <  02}, 


r2  =  {(ai,i/)|6i<t/<62},  (2) 


Ta  =  {(a:,^)  \  ai  <  x  <  02},  r4  =  {(02,2/)  \  h  <  y  <  62}. 

For  simplicity,  we  assume  that  the  controls  on  80,2  are  separable,  i.e.,  they  are  the  product  of  a  function 
of  time  and  a  function  of  the  spatial  variables.  With  this  assumption,  we  specify  conditions  on  80.2  the 
form 

Cl  ^2^ 

w(t,  Ti)  =  ^ui,j(<)<Ai,i(x),  Ti)  =  ^U2,i(t)<^2,i(l/),  (3) 

t=i  «=i 


C3 

U,(t,r3)  =  ^U3,i(«)(^3,i(x),  »x(«,r4)  =  ^U4,i(t)<A4,i(y)- 

i=l  i=l 

In  (3),  Cj  is  the  number  of  controls  on  Tj.  The  control  uj^i  is  the  z-th  control  on  Fj.  The  function  )  is 
a  function  describing  the  influence  of  the  i-th  control  on  Tj. 

On  80i^  we  specify  Dirichlet  and  Neumann  conditions  of  the  form 

w{t,  a,  y)  =  f{y), 


w{t,x,c)  =  0,  w{t,x,d)  =  0. 

In  (4),  f{y)  is  an  inlet  condition  on  the  left,  and  is  analogous  to  the  inflow  condition  specified  in  many  channel 
flow  problems.  The  Neumann  condition  is  an  exit  condition  similar  to  those  specified  in  fluid  problems  where 
the  flow  leaves  the  domain.  The  remaining  Dirichlet  conditions  are  analogous  to  the  no-slip,  no  penetration 
conditions  specified  in  channel  flow  configurations. 

To  complete  the  model,  we  specify  an  initial  condition  of  the  form 

w{0,  X,  y)  =  wq(x,  y)  G  L^(D).  (5) 


The  Control  Problem 

To  facilitate  discussion  of  the  control  problem,  the  initial  boundary  value  problem  given  by  (1),  (3),  (4), 
and  (5)  is  linearized  about  w  =  0.  The  result  is  written  in  abstract  form.  For  the  sake  of  simplicity,  we 
assume  that  the  inlet  condition  f{y)  satisfies  f{y)  =  0  in  the  following  discussion  of  the  abstract  formulation. 

Denote  the  total  number  of  controls  as  Nc^  i  e.  Nc  =  C\  +  (^2  -f-  Ca  -f  (^4.  Specify  F/v  according  to 
F;v  =  •  c  <  2/  <  d},  the  portion  of  80i  on  which  the  Neumann  outflow  condition  is  specified.  Finally, 

let  the  space  V  be  defined  according  to 

V  =  Hl{0)  (^) 

=  {z;  €  H\0)  :  v  \dQ2=  0  =  ^  and  (T m)  =  0}  .  (7) 
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Define  a(*,  )  to  be  the  symmetric  bilinear  form  a{z,v)  =  We  define  the  operator  Aq  : 

V  V'  according  to  [Aqz]  (v)  =  a{z,v). 

Let  W  =  //2(Q)  n  and  define  A-i  \W  H,  where  H  is  the  state  space,  to  be  the  self-adjoint 

restriction  operator  [i4_i]  !/;(•)  =  Wxx{')  +^yy(*)-  By  defining  ||2(-)l|w  as  the  graph  norm 

IN(*)llw  =  II  [^-l]  ^{')\\h  +  ||2;(-)l|//> 

it  follows  that  the  injections  W  C  V  C  H  =  H'  C  V'  C  W'  a^re  all  continuous  and  dense. 

We  now  lift  the  operator  Aq  :  V  ^  V'  to  an  operator  Ai  \  H  W'.  Formally,  we  perform  an  integration 
by  parts  and  define  A\  :  H  W'  according  to 

[Aiz]{w)=  /  z{x,y)Aw{x,y)dx  (9) 

Jq 

for  all  w  €  W'.  ^ 

Utilizing  the  Dirichlet  map  D  :  as  discussed  in  [19],  we  define  the  operator  B  :R  ^  ^  W' 

according  to 

B  =  -AiD.  (10) 

The  linearized  system  is  then  formulated  in  W'  as 

w{t)  =  [A\]w{t) Bu{t)  (11) 


w{0)  =  Wq.  (12) 

We  are  now  in  a  position  to  present  the  control  strategy  utilized  for  the  system  given  by  (1),  (3),  (4), 
and  (5).  We  construct  a  linear  feedback  control  law  from  the  system  given  in  (11)-(12). 

For  a  >  0,  define  the  cost  functional 


^)  / 

^0 


{{Qw,w)h  4-  {Ru,u)u}e^^^dt, 


(13) 


where  Q  is  a  self-adjoint  state  weight  operator  satisfying  Q  >  0.  i?  is  a  self-adjoint  control  weight  operator 
satisfying  R  >  0.  Discussion  and  theoretical  results  for  control  problems  with  the  cost  functional  given  in 
(13)  can  be  found  in  [7,  8]. 

In  this  work,  the  operator  Q  in  (13)  is  specified  to  place  a  larger  weight  over  a  particular  region  of  D. 
Define  Q.q  C  Q  to  be  the  region  in  which  control  of  the  state  is  most  important.  Specify  ^(-^  •)  •  D  — >  R 
according  to 

= 

where  gi,  q2  are  positive  constants.  By  specifying  qi  q2  and  defining  Q  in  accordance  with  (14),  states  in 
Qq  are  weighted  more  heavily  in  the  cost  functional. 

In  a  similar  fashion,  the  control  weight  operator  R  allows  for  larger  weights  on  particular  controls.  Define 


{x,y)  G  n  \  fig, 
{x,y)  e  Dq, 


(14) 


r{u)  according  to 


r{u) 


Vi,  u  a  control  on  Fi, 
r2,  u  a  control  on  r2, 
r3,  u  a  control  on  Fa, 
r4,  u  a  control  on  F4, 


(15) 


where  ri,r2,r3,  and  r^  are  positive  constants.  By  defining  R  in  accordance  with  (15),  controls  on  different 
regions  of  dQ.2  can  be  weighted  differently  in  the  cost  functional. 

Now  that  Q  and  R  are  specified,  the  optimal  control  problem  we  consider  is  to  minimize  (13)  subject  to 

the  constraint 


-—IV  =  Aw  T  Bu, 
dt 


(16) 


n;(0)  =  Wo, 


(17) 
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where  A  and  B  are  the  operators  from  (11).  As  can  be  seen,  we  find  the  feedback  control  law  from  the 
linearization  of  the  nonlinear  problem. 

As  shown  in  [7,  8],  for  an  a-controllable  system,  the  a-LQR  problem  has  a  unique  solution  of  the  form 

Uopt  =  -Kw  (18) 


=  -R-^B*Pw,  (19) 

where  P  is  the  unique,  symmetric,  non-negative  solution  of  the  algebraic  Ricatti  equation 

[A  -f  ayP  -h  P{A  -ha)  -  PBR-^B*P  +  Q  =  0.  (20) 

It  can  be  shown  that  the  action  of  the  gain  operator  K  in  (18)  has  an  integral  representation  of  the  form 

(21) 


Kw{t,x,y)  =  /  k{x,y)w{t,x,y)dxdy, 
JQ 


where  the  functional  gain  /c(-,  •)  G  L^{^)  [9].  The  functional  gains  indicate  how  much  information  each  state 
contributes  to  the  control  and  provide  information  about  optimal  sensor  placement. 

Once  the  gain  operator  K  is  obtained,  the  linear  control  law  is  placed  into  the  nonlinear  system  according 
to 


^w  =  (A-  BK)w  +  G{w)  +  F, 
at 


(22) 


w{0)  =  Wo,  (2^) 

where  G{w)  is  a  nonlinear  operator  resulting  from  the  convective  terms  in  (1)  and  F  is  a  forcing  term 
resulting  from  the  inlet  condition  /(y). 


Numerical  Implementation 


We  utilize  a  finite  difference  approximation  to  numerically  solve  the  initial  boundary  value  problem  with 
and  without  control.  In  the  work  presented  here,  we  specify  a  uniform  grid  on  Cl  with  step-size  h. 

For  the  second  order  derivatives  in  the  x  and  y  directions,  we  implement  the  standard  second  order 
approximations 


/i2 


(24) 


a2 


(a^t.yj)  -  ^  (w'tj+l  -  ■ 


(25) 


Extra  care  is  needed  when  discretizing  tlie  nonlinear  terms.  At  high  Re,  the  convective  terms  dominate 
the  dynamics  of  the  system  and  spurious  oscillations  around  the  shock  occur  if  an  inappropriate  numerical 
scheme  is  cho.sen.  As  discussed  in  [14],  we  utilize  a  mixture  of  central  differences  and  a  donor  cell  discretization 
for  the  convective  terms.  These  approximations  are  of  the  form 


dx 


(26) 


—  W^{Xi,yj)  ~  +«^t,j  +  i)^  -  (w'i.j-l  +  W'i.j)  ] 


+  ^  +  W«,>+l|  (^ij  -  +  ^i,i\  > 


(27) 
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where  7  G  [0, 1]. 

After  applying  the  finite  difference  approximations  given  by  (24)-(27)  and  incorporating  the  conditions 
on  dfti  and  dQ.2^  we  obtain  a  semi-discrete  approximation  of  the  form 

-w^  =  A^w^  +  B‘^u  +  G^  (28) 

dt  ^  ^ 


(29) 

where  we  have  utilized  superscript  N  to  indicate  that  we  have  a  finite-dimensional  approximation  to  the 
infinite-dimensional  distributed  parameter  system. 

The  finite-dimensional  a-LQR  problem  becomes 

w'^fQ{w^)+u'^Ru\e‘^‘^*dt  (30) 


"i"!  [( 


subject  to  ~  +  B^u,  (31) 

«,^(0)  =  Wo^.  (32) 

In  the  finite-dimensional  case,  Q  is  a  diagonal,  symmetric,  positive  semi-definite  matrix  consisting  of  state 
weights,  fl  is  a  diagonal,  symmetric,  positive  definite  matrix  of  control  weights. 

Solving  the  finite-dimensional  a-LQR  problem  yields  a  gain  matrix  ,  resulting  in  an  optimal  control 
of  the  form  Uopt  -  -K^w^ .  The  control  is  placed  into  the  finite-dimensional  system  according  to 

-w^  ^(A^ (33) 
dt  ^  ^  ^ 

w^{0)  =  Wq. 

The  system  (33)-(34)  is  solved  via  a  4th  order  Runge-Kutta  method. 

Numerical  Examples 


Figure  1:  The  inlet  condition  f{y). 

We  now  illustrate  the  effectiveness  of  the  control  strategy.  First,  we  illustrate  the  system  behavior  without 
control  by  specifying  that  the  obstacle  boundary,  80.2 be  held  at  zero  as  time  evolves.  We  then  implement 
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the  control  strategy,  and  present  results,  for  different  values  of  a  in  the  cost  functional.  We  conclude  this 
section  with  a  discussion  of  the  influence  of  the  outflow  boundary  condition  on  the  functional  gains. 

In  all  examples  presented,  the  regions  and  1^2  are  specified  as  =  (0,1.5)  x  (0,  .48)  and  Q2  = 
[.16,  .24]  X  [.12,  .32].  The  step-size  used  in  the  construction  of  the  uniform  grid  is  h  =  .02.  The  inlet 
condition,  /(y),  is  defined  as  in  Figure  1. 

The  values  we  specify  for  K\  and  K2  in  the  governing  equation  are  Ki  =  I  and  K2  =  0.  By  specifying 
Ki  and  K2  in  this  way,  the  inlet  condition  will  result  in  solutions  which  propagate  from  left  to  right.  The 
value  of  7  that  we  specify  in  the  discretization  of  the  nonlinear  terms  is  7  =  1.  In  all  results  presented,  the 
initial  condition  is  defined  by  wo{x,y)  =  0. 

Example  1 

For  the  first  problem,  we  illustrate  the  behavior  of  the  system  in  the  absence  of  control.  This  will  provide 
a  baseline  against  which  to  judge  the  effectiveness  of  the  control  in  the  later  examples.  In  this  example,  we 
set  Re  =  500  and  simulate  the  solution  to  time  T  —  20. 


Figure  2:  The  solution  in  the  absence  of  control. 

As  can  be  seen  in  Figure  2,  the  specification  of  the  inlet  condition  f{y)  results  in  a  solution  propagating 
from  left  to  right.  In  particular,  the  solution  is  a  traveling  shock.  The  shock  encounters  the  obstacle,  and 
a  portion  of  it  propagates  above  and  below.  As  time  evolves,  solution  values  increase  downstream  of  the 
obstacle  until  a  steady-state  solution  is  reached. 

Example  2 

We  now  implement  boundary  control  on  0.2-  In  particular,  we  specify  that  there  are  two  controls  each  on 
the  front  and  back  of  the  obstacle.  Similarly,  we  specify  that  there  is  one  control  each  on  the  bottom  and  top 
of  the  obstacle.  For  simplicity,  we  specify  that  the  control  influence  functions  {0j,i(  )}  are  piecewise  constant. 
The  control  weights  ri,r2,r’3,  and  in  (15)  are  specified  to  be  500,  100,  500,  and  1000,  respectively. 

The  control  objective  is  to  drive  nonzero  solution  values  in  a  region  downstream  of  the  obstacle  to  zero. 
To  this  end,  we  specify 

Qq  =  [.5,  .75]  X  [.16,  .32].  (35) 
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To  severely  penalize  nonzero  states  in  Qg,  we  specify  that  q\  =  I  and  q2  =  1000.  These  values  result  in  states 
in  fig  being  weighted  much  more  heavily  in  the  cost  functional  u{t))-  In  order  to  make  a  comparison 

of  the  impact  of  a  on  the  effectiveness  of  the  control,  we  specify  that  a  =  0  in  this  example.  Of  course, 
a  =  0  corresponds  to  the  traditional  LQR  control  formulation. 


Figure  3:  The  steady-state  solution  for  a  =  0. 


As  can  be  seen  in  Figure  3,  the  traditional  LQR  formulation  of  the  control  problem  is  not  very  effective 
in  this  case.  Nonzero  values  in  Qg  are  only  slightly  reduced,  even  though  a  severe  weight  is  placed  on  these 
values  in  the  cost  functional.  The  heavy  penalty  on  states  in  Q.q  has  a  direct  influence  on  the  functional 
gains,  as  seen  in  Figure  4.  The  functional  gains  for  controls  on  the  bottom,  top,  and  rear  of  the  obstacle  have 
large  magnitudes  over  Qq,  As  a  result,  states  in  f)g  contribute  significantly  to  those  controls.  Conversely, 
states  in  Qg  do  not  have  a  significant  contribution  to  controls  on  the  front  of  the  obstacle.  For  controls  on 
the  front,  the  functional  gains  have  large  magnitudes  near  the  control  location. 
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Figure  4:  Results  for  q  =  0.  Functional  Gains  (left),  Controls  (right). 
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Controls  resulting  from  the  traditional  LQR  formulation  have  very  small  magnitude,  as  evident  in  Figure 
4.  They  are  not  effective  at  driving  nonzero  states  in  Hq  to  zero. 

Example  3 

We  now  investigate  the  effectiveness  of  the  control  in  the  case  of  nonzero  a.  In  this  example,  the 
parameters  defining  the  problem  are  specified  as  in  Example  2  with  the  exception  that  we  set  a  =  0.2.  As 


Figure  5:  The  steady-state  solution  for  a  =  0.2. 


can  be  seen  in  Figure  5,  the  control  is  much  more  effective  at  suppressing  nonzero  states  in  Hq  for  the  case 
a  =  0.2.  A  majority  of  the  nonzero  states  in  D.q  are  greatly  reduced  when  compared  to  their  corresponding 
values  in  Figure  2.  The  control  is  much  more  effective  in  this  case. 


Figure  6:  Results  for  a  =  0.2.  Functional  Gains  (left).  Controls  (right). 
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Specifying  nonzero  a  has  a  direct  impact  on  the  functional  gains.  From  Figure  6,  the  functional  gains 
have  much  larger  magnitudes  than  was  the  case  for  q  =  0.  In  addition,  the  gains  are  large  for  states  in 
Hq,  as  well  as  for  many  states  outside  of  fig.  As  a  result,  a  majority  of  states  downstream  of  the  obstacle 
contribute  to  the  control.  The  resulting  controls  are  of  larger  magnitude  than  was  the  case  for  a  =  0,  as 
evident  in  Figure  6.  However,  the  magnitudes  of  the  controls  are  not  unreasonably  large.  There  is  a  tradeoff 
in  the  construction  of  Q,  /?,  and  the  value  specified  for  a.  Obviously,  the  state  and  control  weights,  as  well 
as  the  parameter  a,  need  to  be  chosen  carefully  in  order  to  yield  a  control  that  is  effective,  but  not  too  large. 

Curiously,  all  of  the  functional  gains  of  this  example  have  their  maximum  values  at  the  outflow  on  the 
right  of  the  domain.  At  the  right  boundary,  a  Neumann  condition  is  specified  in  an  analogous  way  to 
the  outflow  condition  specified  in  many  flow  problems  where  the  flow  leaves  the  domain.  This  condition 
is  specified  in  order  to  allow  the  numerical  solution  to  freely  pass  through  the  right  boundary.  From  this 
example,  it  appears  that  a  large  portion  of  the  control  effort  is  directed  at  forcing  the  numerical  outflow 
condition.  We  further  investigate  this  issue  in  the  next  example. 

Example  4 

For  the  final  example,  we  investigate  the  influence  of  the  outflow  boundary  condition  on  the  functional 
gains.  As  seen  in  the  previous  example,  it  appears  that  much  of  the  control  effort  is  aimed  at  forcing  the 
Neumann  condition  specified  at  the  right  boundary.  To  illustrate  the  influence  of  the  outflow  condition  on 
the  functional  gains,  in  this  example  we  specify  all  parameters  as  they  were  in  Example  3  with  the  exception 
that  a  Dirichlet  condition  is  specified  on  the  right  boundary  of  the  form 

ti;(<,rN)  =  0.  (36) 
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Figure  7:  Results  for  a  =  0.2,  Dirichlet  Condition.  Functional  Gains  (left),  Controls  (right). 


In  the  case  of  a  =  0,  modifying  the  boundary  condition  along  has  no  effect  on  the  functional  gains. 
The  gains  obtained  with  the  condition  specified  in  (36)  are  identical  to  those  obtained  in  the  Neumann 
case,  i.e.,  those  of  Figure  4.  However,  modifying  the  right  boundary  condition  has  a  significant  effect  on  the 
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functional  gains  for  the  case  a  =  0.2,  as  evident  by  comparing  the  gains  in  Figure  7  with  those  of  Figure 
6.  The  gains  computed  with  the  Dirichlet  condition  given  in  (36)  have  much  smaller  magnitude  and  obtain 
their  maxima  on  the  interior  of  the  domain.  Moreover,  the  controls  obtained  in  the  Dirichlet  case  have 
much  smaller  amplitudes  than  they  did  in  the  Neumann  case.  Clearly,  the  outflow  condition  needs  to  be 
chosen  carefully.  The  outflow  condition  is  primarily  specified  as  a  numerical  convenience  in  order  to  allow 
the  computed  solution  to  freely  exit  the  computational  domain.  However,  as  evident  in  this  example,  the 
condition  specified  for  outflow  can  have  a  dramatic  eflect  on  the  functional  gains  and  the  resulting  controls, 
depending  on  the  cost  functional  used  in  the  formulation  of  the  control  problem. 


TTm  Steady  State  Sotuiion 


Figure  8:  Dirichlet  gains  applied  with  Neumann  outflow  for  a  =  0.2. 

To  verify  that  a  large  portion  of  the  control  effort  of  Example  3  is  aimed  at  forcing  the  Neumann  outflow 
condition,  we  apply  the  functional  gains  of  this  example  to  the  system  of  Example  3  where  a  Neumann 
outflow  condition  is  specified.  The  resulting  steady-state  solution  is  shown  is  Figure  8.  The  steady-state 
solutions  shown  in  Figure  5  and  Figure  8  are  practically  identical.  The  gains  obtained  in  this  example, 
and  their  corresponding  controls,  are  very  effective  at  suppressing  nonzero  states  in  Dq  for  the  system  of 
Example  3.  Moreover,  the  functional  gains  and  controls  of  this  example  have  much  smaller  magnitudes  than 
their  counterparts  in  Example  3.  The  controls  of  this  example  attain  the  control  objective  with  much  less 
effort.  Clearly,  careful  consideration  of  the  outflow  condition  must  be  done  in  order  to  avoid  controls  that 
are  unnecessarily  large  for  the  sake  of  a  convenient  numerical  outflow  condition. 

Conclusions  and  Future  Work 

In  this  paper,  we  have  investigated  the  effectiveness  of  a  particular  linear  feedback  control  strategy 
when  applied  to  a  system  where  the  dynamics  are  governed  by  the  two-dimensional  Burgers’  equation.  We 
specifically  constructed  the  system  in  such  a  way  as  to  incorporate  analogous  conditions  as  those  specified 
in  many  fluid  flow  configurations.  The  strategy  chosen  to  control  the  system  was  very  effective.  However, 
as  seen  in  the  results,  care  must  be  taken  in  specifying  quantities  used  in  the  control  formulation.  In 
particular,  the  construction  of  the  state  and  control  weights  Q  and  R  must  be  done  carefully  in  order  to 
place  appropriate  emphasis  on  regions  of  interest,  and  to  correctly  formulate  the  control  objective.  Moreover, 
the  value  specified  for  the  additional  performance  p)arameter  a  needs  to  be  chosen  carefully  to  ensure  that 
the  resulting  control  is  not  prohibitively  large  for  real-world  applications. 

As  seen  in  the  computational  results,  the  numerical  outflow  condition  can  have  a  dramatic  effect  on  the 
functional  gains  and  their  corresponding  controls.  The  outflow  condition  needs  to  be  chosen  carefully  so  that 
unnecessary  control  effort  is  not  expended  solely  for  the  sake  of  numerical  convenience. 

As  our  motivation  for  considering  the  problem  presented  in  this  paper  was  an  investigation  into  the  utility 
of  linear  control  techniques  to  fluid  flow  problems,  future  work  includes  investigating  these  techniques  for  the 
Navier-Stokes  equations.  There  are  several  issues  that  must  be  addressed  in  that  situation.  The  case  of  large 
Reynolds  number  in  the  Navier-Stokes  equations  leads  to  systems  where  the  dynamics  are  dominated  by 
nonlinearities.  More  energy  is  required  to  obtain  the  control  objective  if  the  control  enters  into  the  problem 
linearly.  One  aspect  of  future  work  involves  developing  nonlinear  feedback  control  methods  to  overcome  this 
difficulty. 
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Obviously,  the  large  number  of  grid  points  needed  to  resolve  a  flow  field  is  a  significant  obstacle  in 
implementing  the  control  techniques  in  this  paper  to  a  fluid  flow  problem.  In  this  paper,  we  investigated 
the  effectiveness  of  the  control  under  optimal  circumstances,  .i.e.,  full  state  information  was  available  for 
feedback.  However,  full  state  information  will  not  be  made  available  in  a  fluid  flow  problem  as  the  resulting 
control  problem  will  be  prohibitively  large.  As  a  result,  an  investigation  into  optimal  sensor  placement  is 
also  an  area  of  future  work. 

Finally,  more  sophisticated  methods  of  finding  the  gain  operator  should  be  utilized  in  the  case  of  optimal 
control  of  fluid  flow.  In  that  case,  finding  the  gains  via  an  algebraic  Ricatti  equation  is  prohibitively  expensive. 
Methods  exploiting  the  structure  of  the  problem,  e.g.  methods  based  on  the  Chandrasekhar  equations  [6,  17], 
have  the  potential  of  providing  significant  savings  in  the  case  of  flow  control. 
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